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«he  bulk  crystal.  The  analysis  of  the  Bragg-peak  intensities  is  carried  out 
based  on  the  usual  kinematical  theory  and  the  dynamical  (two-beam  approxi- 
mation) theory  by  Blackman.  TtW  fonsar  theory  fails  completely  in  accounting 
for  small  angle  reflections  (s(y5(j^),  whereas  the  latter  explains 
successfully  the  whole  observed  intensities,  except  for  the  2nd  and  3rd 
order  reflections  (400),  (440),  (333)  and  (600).  From  the  damping  of  peak 
intensitiesj^he  temperatures  of  these  clusters  are  estimated  to  be  of  the 
order  of  iMW.  » ^ 


A close  ins^e\tion  on  the  comparison  between  the  observed  and  theoreti-  . 
cal  intensities,  howver,  reveals  periodic  fluctuations  in  the  observed  peak 
heights  as  a functio^of  the  scattering  angle,  indicating  a superposition 
of  some  diffuse  scatt»ing.  The  amplitudes  of  the  fluctuations  are  more 
marked  for  the  smaller  size  clusters.  The  characteristics  of  the  fluctuation! 
are  accounted  for  by  introducing  llquld-like  random  atomic  configurations 
into  the  microcrystalline.  This  model  explains  also  the  abnormal  Bragg-peak 
intensities  of  Ar  clusters,  reported  in  the  literature. 
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INTRODUCTION 

Little  has  been  known  about  the  structure  of  micro crystals, 
despite  extensive  theoretical  and  experimental  works.  Electron 
diffraction  observations  for  clusters  above  a certain  size  show 
the  Bragg  reflection  peaks,  chsuract eristic  of  the  regular  lattice  of 
bulk  crystals.  Does  this  mean  that  such  a cluster  is  simply  a micro- 
replica of  bulk  materials?  The  question  could  be  answered  by  a detailed 
analysis  of  the  reflection  intensity,  instead  of  the  simple  lattice 
Parameter  analysis*  Farges  et  al.^^  observed  the  anomaly  in  some 
of  the  reflection  intensities  of  Ar  clusters  consisting  of  500'*^1000 
atoms,  suggesting  an  amorphous  nature  of  the  micrdcrystal  structure. 

For  very  small  sizes  of  clusters  tsay,  <C  500  atoms  ),  diffraction 
patterns  become  similar  to  those  in  bulk  liquid,  as  observed  in  Ar 
clusters.^' According  to  theoretical  model  calculations, the 
structure  of  microclusters  is  different  in  quality  from  the  bulk  lattice 
structure  and  even  amorphous  for  small  sizes  of  aggregates.  The 
structure  and  stability  of  microclusters  depend  strongly  on  the  size 
of  aggregates.  Microclusters  are  in  a pseudo-stc^dy  state  and  their 
structures  evolve  through  a so-called  critical  size.  Such  a dynamical 
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aspect  of  the  cluster  structure  is  an  important  clue  to  understand 

phase  change  phenomena,  nucleation  kinetics,  and  catalytic  activities 

1) 

of  metal  clusters,  etc.  ' 

Thus,  our  immediate  concern  is  to  shed  some  light  on  such  a cluster 
structure  hy  means  of  molecular  beam  - electron  diffraction.  So  far 

g) 

metal  clusters  have  not  been  studied  by  this  method,  though  thin 
film  experiments  are  well  known.  An  advantage  of  this  technique  over 
the  thin  film  is  that  clusters  can  be  produced  by  homogeneous  nucleation, 
and  their  diffraction  patterns  are  free  from  the  sample-supporting 
material  for  films.  Such  supporting  materials,  being  amorphous, 
give  a diffuse  scattering  background,  by  which  the  measurement  of 

reflection  peak  intensities  is  sometimes  disturbed.  A difficulty  in 

• 

the  present  method,  however,  is  to  prepare  the  cluster  sample  with 
desired  sizes.  Sv^ersonic  expansions,  as  successfully  used  for  simple 
gas  samples  (Ar,  COg,  etc.),’  * ’ are  not  practical  for  metals. 

Cooling  of  metal  vapor  by  a cold  carrier  gas  produces  easily  metal 
clusters,  but  for  the  sampling  of  desired  sizes  of  clusters  a proper 
nozzle  must  be  designed  f>r  the  diffraction  experiment.  A simple 
nozzle  for  this  purpose  has  been  constructed  to  produce  a metal  cluster 
beam*^  This  nozzle  has  enabled  us  to  study  the  structure  of  lead  metal 
clusters,  Itr  the  first  time,  in  their  free  state  (gas  phase). 

Lead  has  been  chosen  ^ the  following  reasons i 1)  technical 
ease  of  vaporization i 2)  well-known  bulk  properties  which  are 
useful  for  the  analysis i 3)  a simple  structture  of  the  bulk  crystal 
(face -centered  cubici  f.c*o«),  which  is  the  same  as  Ar  solid)——  Ar 
clusters  are  most  extensively  studied  and  standard  samples  for  the 
cluster-structure  studyi  and  4)  possible  importance  of  dynamical 
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X 10) 

(multiple  scattering)  effects  due  to  the  heavy  atoms  (Z  * 82),  being 
different  from  the  case  of  Ar  clusters. 


t . 


EXPERIMENTAL 

The  metal-cluster-beam  nozzle  used  in  this  study  is  schematically 
shown  in  Pig.  1|  a more  detailed  description  of  this  nozzle  will  be 
reported  elsewhere.®^  Lead  (99.999?^  purity,  pvirchased  from  Goldsmith 
Chem.  & Metal  Corp.)  was  evaporated  from  an  electrically  heated  tungsten 
filament  oven  (H) , and  mixed  with  argon  carrier  gas  at  about  0.6  cm 
downstream  from  H.  The  ambient  temperature  (Tg^)  around  this  mixing 
region  was  monitored  with  a Chromel-Alvimel  thermocouple  (TC) , while 
the  temperature  of  metal  on  the  filament  was  measxired  by  an  optical 

pyrometer  through  a small  back  window  (HP).  No  special  cooling  system 
was  required  for  this  nozzle  assembly  because  of  a small  heat  capacity 
of  the  filament  oven,  though  the  filament  temperature  was  typically 
900 —lOOO^C*  Condensation  of  metal  vapor  occurs  through  an  isobaric 
cooling  by  surgon  carrier  gas  at  the  mixing  region,  where  the  distance 
(I^)  from  the  heater  H.  to  the  first  nozzle  aperture  (N^)  is  adjustable 
(1.6  — 5.6  cm).  The  condensates  further  cool  by  two-stage  adiabatic 
expansions  throu^  (diami  0.076  cm)  and  Ng  (diami  0.050  cm)  nozzle 
apertures.  The  pressure  drop  by  these  free  jet  expansions  was  of  the 
order  of  2zl0"5,  The  expansions  form  a well- collimated  metal  cluster 
beam  with  an  angular  width  of  about  5°i  at  about  50  cm  downstream,  a 
circular  shape  deposit  (c.a.  4.5  cm  in  diameter)  of  metal  was  observed. 
In  order  to  change  an  average  size  of  metal  clusters,  the  distance  L^, 
argon  gas  pressure  (P^)  and  the  filament  temperature  (Tg^)  were 
varied.  Although  no.  attempt  was  made  to  study  the  nucleatlon  kinetics, 


it  was  found  by  some  trial- and- error  experiments  that  Pa  and  were 
effective  parameters  to  control  the  clxister  size.^^ 

The  metal  cluster  beam  generated  by  this  nozzle  was  crossed  with 
an  electron  beam^®^(with  an  accelerating  voltage  of  40  kV  and  an 
electron  beam  current  of  4,6  jtA)  at  0.o64  cm  distant  from  the  nozzle 
tip  (N2)  in  the  diffraction  chamber,  where  the  ambient  pressure  of 
(0.8  — 1,2)x10"5  torr  was  maintained  during  the  experiments.  The 
metal  beam  diameter  at  the  point  to  cross  the  electron  beam  was  measured 
to  be  0,056  cm  from  a size  of  a trace  of  metal  deposit  on  a thin  A1 
plate  which  was  placed  at  the  crossing  point.  This  spread  of  beams 
causes  only  a minor  correction  to  the  line  broadening  of  diffraction 
patterns.  Diffraction  patterns  of  Pb  cluster  beams  were  recorded  on 
Kodak  Electron  Image  plates  (5-^X4")  at  a nozzle-to-plate  distance  (L) 
of  53  cm  without  a sector.  The  magnetic  field  from  the  heater  current 
(D,C,  regulatedi  max,  30Amp)  produced  on  the  photographic  plate  a 
deflection  of  less  than  0,2  cm,  which  was  compensated  by  moving  the 
beam  stop  and/or  tilting  the  electron  gun  unit  mechanically.  For  the 
instrumental  calibration,  diffraction  patterns  for  standard  thin  metal 
films  of  Al,  Au  and  Pb  were  also  recorded.  These  provided  the  camera 
constant  LX  (X  ■ electron  wave  length)  and  an  estimate  of  the  instrumental 
line  broadening  of  Debye-Scherrer  rings. 

Three  typical  plates  (denoted  as  A,  B and  C)  were  chosen  for  the 
following  analyeis.  The  experimental  conditions  are  given  in  Table  1, 

The  densitometric  traces  for  these  plates  are  shown  in  Pig,  2,  In 
order  to  convert  the  photographic  density  (d)  to  the  relative  intensity 
(I),  the  following  equation  was  used. 
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I -•  d(l+  cd)  (1) 

, where  c « 0.1  for  d^2.0.  This  equation  fitted  a measured  density- 
exposure  calibration  curve. 


ANALYSIS  AND  RESULTS 
I)  Peak  Positions  and  Widths 

According  to  the  Bragg  condition,  the  peak  position  of  diffraction 
pattern  * ven  by, 

8 a 2Tt/  (2) 

.where  s * (47t/>.)sin9  (diffraction  angle  *29),  and  is  the  inter- 

planar  distance  with  the  Miller  indices  (hkl) . For  a f.c.c.  crystal 
and  a small  diffraction  angle  (in  the  present  case,  9 < 4. 5*10 "^radians ) , 
eq.  (2)  can  be  expressed  by  a radius  (r)  of  the  Debye-Scherrer  ring 
for  plane  (hJcl) , the  camera  constant  (LA)  and  the  lattice  constant  (a). 

r/LXot  1/dj^^  » (h2  4-  k^  ♦ (3) 

Once  each  ring  is  indexed  by  the  standard  procedure  (see  Pig.  2), 

the  measuronent  of  the  ring  diameter  provides  the  lattice  parameter. 

The  systematic  error  due  to  the  small  angle  approximation  in  eq.  (3) 

can  be  sufficiently  corrected  by  the  use  of  r - (L  * 53  cm) 

instead  of  r.  This  correction  becomes  impoirtant  for  outer  ringsi  e.g. , | 

0 i 

the  correction  for  a (531)  ring  results  in  an  increase  of  0.01  A in 
the  lattice  parameter.  The  results  are  sxanmarized  in  Table  2 for  the 
three  chosen  plates  A,  B and  C. 

t- 
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The  width  of  Bragg  peaks  provides  an  estimate  of  the  average 
cluster  size,  as  has  been  used  in  X-ray  powder  patterns}^^  after 
correcting  the  instrumental  broadening,  which  was  estimated  by  taking 
diffraction  patterns  of  standard  thin  films i see  experimental  section. 

The  broadening  due  to  the  cluster  beam  spread  was  of  negligible  im- 
portance, compared  to  the  above  source.  For  the  estimate  of  the 

average  cluster  size  (diameten  D) , weUrresolved  (111)  and  (220)  peaks 

12) 

were  used  with  assuming  the  unit  shape  factor,  and  the  results  are 
given  in  the  lower  section  of  Table  2 for  the  three  plates.  The  errors 
given  in  the  Table  were  estimated  from  uncertainties  in  the  instrumental- 
broadening correction  and  in  the  line-width  measurement.  The  average 
size  D can  also  be  obtained  independently  from  the  following  peak- 
intensity  analysis. 

II)  Peak  Intensity 

The  relative  intensity  of  Bragg  peaks  was  obtained  by  subtracting 
the  background  intensity  from  the  experimental  curves,  which  were 
calibrated  by  eq.  (1)  to  the  relative  intensities.  Smooth  background 
curves  ( base  line  for  the  Bragg  peaks  ) were  drawn  without  ambiguity 
for  the  three  plates.  In  order  to  examine  systematic  background  errors, 
if  any,  due  to  the  instrument,  diffraction  patterns  of  pure  argon  gas 
were  taken  under  experimental  conditions  similar  to  those  in  the  metal 
sample.  No  discernible  systematic  fluctuatioa  or  periodicity 
was  observed.  The  observed  peak  intensities  lots*  normalized 

to  100  at  the  (111)  peak,  are  listed  in  Table  2. 

The  relative  peak  intensity,  Isyagg*  calculated  by  the 

usual  kinsmatical  expression  for  X-ray  powder  pattem^^^  with  a small 
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angle  approximation! 

^Bragg  “ |''h]al^‘’’hW.-^?lcl-e-2"  ^ (4) 

l^ial  =Z'J  . exp[2iti(hxj  + ky^  f IZj)]  = 4-f  (5) 

(Where  fj  is  the  elastic  atomic  scattering  factor  of  electron  for 
the  atom  with  a coordinate  ^ unit  cell.  is 

a multiplicity  of  the  (hkl)  reflection  with  the  interplanar  distance 
'^hkl*  Debye-Waller  factor  e”^  can  be  expressed  with  an  isotropic 

square  mean  amplitude  iv^)  of  lattice  vibrations!  = exp(-n^s^). 

A plot  of  l>^tiobs/^kin^  should  give  a straight  line  v/ith  a 

negative  slope  of  u^,  if  eq.  (4)  holds  for  the  present  case.  However, 
such  plots  showed  a marked  deviation  from  a straight  line  for  low 
angle  peaks  ( s ^5  and  for  the  2nd  order  reflections  (4oo)  and 

(440).  The  observed  intensities  for  these  peaks  are  too  weaki  e.g., 

80?4  and  50^  lower  than  the  kinematical  values  for  the  (ill)  and  (400) 
peaks,  respectively.  This  is  illustrated  in  Fig.  3 by  the  curve  (5) 
in  the  case  of  Plate  A.  The  failure  of  the  kinematical  theory  in 
electron  diffraction  is  to  be  expected  for  a micro crystal,  specially  one 
which  consists  of  heavy  atoms  such  as  Pb  (Z  » 82) , since  the  attenuation 
(or  extinction)  of  the  original  electron  wave  due  to  interatomic 
multiple  scatterings  (or  dynamical  interaction  effects)  is  more  im- 
portant than  that  of  X-ray,  because  of  the  stronger  interaction  between 
electrons  and  atoms. 

For  the  correction  of  intensities  due  to  the  dynamical  effects, 

14)  • 

the  Blackman  formula,  known  as  the  two-beam  (the  incident  and  one 
diffracted  beam)  approximation,  is  commonly  used  in  thin  film 
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experiments  by  electron  diffraction^’^ • 17)  it  is  kno\vn  that  this 

theory  is  only  applicable  fororystzis  containing  one  kind  of  atom  or 
else  containing  only  light  atoms,  and  that  it  is  inadequate  for  the 
higher  order  reflections. 

According  to  the  Blackman  theory,  the  dynamical  peak  intensity 
Udyn)f  corresponding  to  in  eq.  (4),  is  given  by 

^dyn  * (6) 

A = (2mtig/h^)*  jPhkll 

.where  Jq(2x)  is  the  zero  order  Bessel  function,  j^hkll  kinematicai 

structure  factor  given  in  eq.  (5) , D is  the  thickness  of crystals  Cor 

the  diameter  of  a cluster),  and  A is  the  electron  wave  length  tin  the 

. 0, 

present  case.X^  0*061  A). 

Plots  of  ln(lQ^g/I(iyn)  vs,  s^  were  carried  out  with  varying  the 

size  parameter  D systematically,  since  eq.  (6)  depends  on  D;  typical 

plots  for  Plate  A are  shown  in  Fig.  5*<ienoted  as  (1)  — (4).  The 

deviation  of  the  low  angle  reflections  from  a straight  line,  observed 

in  the  previous  kinematicai  analysis,  was  completely  corrected  in  the 

case  of  D = 80  ± 5 ^ with  the  use  of  partial  -wave  elastic  scattering 

factors^®^  for  f in  eq.  (5).  The  size  parameter  thus  obtained  is  in 

excellent  agreement  with  82  ± 10  A,  estimated  from  the  line  widths  ( 

see  the  previous  section).  A similar  analysis  using  the  Born  atomic 

18 ) 

scattering  factors  provided  also  a straight  line,  but  led  to  a too 
small  value  in  D of  6o±3  X (see  (5)  in  Fig.  3),  indicating  the  im- 
portance of  correction  for  intra-etomic  multiple  scattering  effects, 
which  are  taken  into  account  in  the  partial  wave  scattering  factors. 
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As  stated  earlier,  the  2nd  order  reflections,  (400)  and  (440) 

except  (222),  can  not  be  explained  by  the  two-beam  approximation; 

see  Pig. 5.  The  Jrd  order  reflections  (353)  and  (600) , which  coincide 

with  the  first  order  (511)  and  (442)  reflections,  respectively,  were 

also  too  weak.  In  Fig,  3,  the  3rd  order  reflections  were  ignored; 

when  they  are  included,  all  points  for  (511)  and  (442)  in  the  Figure 

should  be  shifted  do\mward  by  0.29  and  0.22,  respectively. 

Parallel  analyses  for  Plates  B and  C were  carried  out  using  the 

partial  wave  atomic  form  factori’®^  Plots  of  Indg^g/Ij^in)  and 

ln(I  V /Ij  ) showed  a behavior  similar  to  the  case  of  Plate  A.  The 
obs'  dyn 

best  straight  lines  were  obtained  with  0 =6o±5  ^ for  Plate  B and 
o 

D = 5o±5  A for  Plate  C,  respectively,  which  are  also  consistent  with 

those  obtained  from  the  line-width  analysis.  These  plots  are  shown 

0 

in  Pig,  4,  compared  with  Plate  A (D  » 80  A) i the  3rd  order  reflections 
are  not  included. 

The  slopes,  u , of  these  plots  were  determined  to  be  0.015  A (Plate 
A),  0.017i  A^  (B)  and  0.012  A^  (C) . Within  the  framework  of  kinematical 
theory  for  lattice  vibration  effects,  can  be  interpreted  as  mean 
(isotropic)  square  displacements  of  atoms,  and  related  to  the  Debye 
characteristic  temperature  0 and  the  sample  temperature  T,  being 
approximated  by  (0  T)*’^’^^^ 

* 5h2T/4  7t^Mk®^  (7) 

.where  h is  the  Planck  constant,  M is  the  atomic  mass  of  Pb,  and  k 

is  the  Boltzmann  constant.  The  Debye  temperature  for  bulk  Pb  crystal  in 

. o 20) 

the  literature  ranges  from  68  to  105  K,  and  for  the  surface  atoms 

it  has  been  reported  to  be  55  ^ 10° K.  By  allowing  a range  of 

0 = 55  — -105°K,  eq.  (7)  gave  a rou|^  estimate  of  the  cluster  tern-  ; 

i 
( 
e 

t 
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peraturei  T (Plate  A)  * 150  * 87°K,  T(B)  = 170  t 90Ok,  and  T(C)  = 

120  t 7o°K. 

Although  the  observed  (first  order)  Bragg  reflections  are  well 

explained  by  the  use  of  the  effective  atomic  scattering  factor  based 

on  the  two-bean  approximation  and  the  kinematical  Debye-Waller  factor, 

a close  examination  of  scattered  experimental  points  in  Fig.  4 shows 

systematic  fluctuations  along  the  straight  lines.  The  amplitude  of 

the  periodic  fluctuations  is  enhanced  as  the  size  of  clusters  becomes 

0 

small  I in  the  case  of  the  largest  size  (Plate  A,  D = OO  A),  such 
fluctuations  die  down  into  the  e9()ezimental  imdora  noise.  This  fact 
implies  that  the  origin  of  the  fluctuations  should  arise  from  size 
effects,  i.e.,  the  nature  of  cluster  structure,  instead  of  other  possible 
sources  such  as  multiple  scattering  effects.  The  observed  fluctuations 
of  Bragg  intensities  are  analyzed  in  the  following  section  using  a 
liquid  model,  and  f\urther  discussed  in  a later  section. 

Ill) 

In  order  to  see  the  fluctuations  clearly,  the  following  function 
was  defined  I 

”<•)  • WlBr.«g 

.Where  I-  * with  the  notations  used  previously. 

Bragg  “y« 

R(s)  for  Plate  C.  in  which  the  most  marked  fluctuations  were  observed, 
is  plotted  in  Pig*5i  the  period  of  R(s)  is  about  s «2  A . This 
slowly  varying  function  of  s suggests  a superposition  of  a typical 
diffuse  scattering  as  observed  in  a liquid  or  gas  target. 

When  the  Bragg  intensity  is  superimposed  on  the  diffuse  scatter- 
ing intensity  Iq,  due  to  the  disordered  atomic  oonf iguratione , the 
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R(s)  in  9q«  (8)  can  be  equated  by  . 

R(s)  ~ 1 ♦ (9) 

.where  k is  a measure  of  the  "disordernessr  but  its  physical  meaning 
is  not  clear  because  of  the  simplicity  of  model  used  herei  in  the 
following  analysis  k is  regarded  merely  as  a normalization  constant. 

Ijj  is  the  Fourier  transform  of  the  atomic  pair  correlation  function 
h(r) , 

Ijj  ~ f^n  Jh(r)exp(i?*f)d?  » f^nH(s)  ( 10) 

.where  f is  the  effective  atomic  scattering  factor  for  Pb  used  in  the 
calculation  of  n is  a mean  number  density,  H(s)  can  be 

expressed  by  the  direct  correlation  function  C(s)  in  momentixm  space, 
after  performing  the  Fourier  transformation  of  the  Or nstein-Zernike 
direct  correlation  function  c(r)|^^ 

H(s)  • C(s)/(1  - nC(s))^  (11) 

C(s)  » jo(r)exp(i?*7)d?. 

In  order  to  evaluate  eq,  (11),  a simple  "liquid  model"  with  particles 
interacting  throu^  the  hard-sphere  potential  was  employed i for  most 
metal  liquids  such  a model  has  been  successfully  applied  to  explain 
the  liquid  structure  factor?^^  According  to  Wertheim^^^  and  Thiele?^ ^ 
the  direct  correlation  function  o(r)  for  the  hard-sphere  potential 
isphers  disaster  6)  is  given  within  the  Percus-Yevick  approximation^^) 
by  a polynoainal* 
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- c(r)  - aL  * (3(r/ff)  ♦ 1(r/<fp 
a « (1  ♦ 2h)2/(1  - 
/5«  -6»iU  4.  ij/2)2/(1  - n)^ 
y-  t^(l  + 2ti)2/2(l  - 

(packing-density  parameter). 

This  gives  an  analytical  form  to  H(s)  in  eq.  (11). 

nH(s)  = [l  ♦ (24tj/x^)  (aj^sinx  - a2Cosx  + - 1 

X SB  <rs 

a3^=a+2p+4r-  24r/x^ 
a^  *s  (A  ♦ ^ + y )x  - 12y/x  + aj 
a,  * -2p/x  * ZAX/y^. 


Then,  R(s)  becomes  the  following  analytical  equation,  using  eqs.  (9), 
(10)  and  (13). 

R(s)  =1+  (iC6^**/<ihki^hkl^{(f  ^ (24ij/x?)  (sj^sinx  - 

agcosx  ♦ aj))’^  - l}^  (l^ 

o * 

Using  O'  « 3.15  A,  2U  k0.012s^  and  = 0.45  (a  typical  value  for 
most  metal  liquids) eq.  (14)  was  calculated  at  the  Bragg  reflection 
angles.  The  results,  being  normalized  to  1.26  (observed  value)  at 
the  (111)  refleotion,  are  plotted  in  Pig.  5,  in  comparison  with  the 
observed  points.  The  observed  periodicity  and  relative  amplitudes 
are  well  reproduced.  Another  choice  of  the  magnitudes  in  O’  and 
could  give  a better  fit,  but  such  an  optimization  was  not  performed. 
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Discussion 

The  lattice  parameter  is  the  moat  direct  information  of  the 
hulk  structure  of  microclusters » whenever  it  is  well  defined,  i.e. 
experimentally^  the  indexing  of  Bragg  reflections  is  possible  (as 
in  the  present  case).  Diffraction  experiments  of  thin  films  indicate 
that  the  lattice  parameter  of  microclusters  is  smaller  than  that  of 
bulk  materials,  and  decreases  with  decreasing  cluster  sizef  ’ ' 

It  is  q,uite  conceivable  that  some  surface  effects  are  operating  in  the 
size  dependence.  One  of  the  plausible  explanations  for  the  observation 
is  a phenomenological  one  to  use  the  concept  of  the  elastic  compression 
by  the  surface  tension  (S).  By  assuming  a spherical  cluster  with  a 
diameter  (D)  and  the  bulk  modulus  (k)/  the  shrinkage  ( ^a)  in  the 
lattice  parameter  (a)  is  approximated  by, 

^a  a -4aS/3kD.  (15) 

When  we  apply  this  equation  to  the  present  case  (Pb) , for  the  estimate 

of  the  order  of  magnitudes,  by  the  use  of  k » 5 XIO^^  dyn/cm^, 

S ai  500  dyn/cm^  and  a » 4.9505  A (of  the  bulk  crystal),  the 

o 

decreased  lattice  parameters  for  D = 80,  6o  and  40  A are  4.942,  4.939 
o 

and  4.934  A,respectively,  being  compatible  with  the  corresponding  observed 
values  I 4.939,  4.938  and  4.935  A (see  Table  2). 

Another  macroscopic  explanation  for  the  present  case  is  due  to  a 
temperature  effect.  The  linear  thermal  expansion  coefficient  for  Pb 
is  27.08 X10“V°C  for  the  range  from  -183°C  to  14°C.  The  lattice 
parameters  for  Plates  A,.  B and  C due  to  this  effect,  using  the  estimated 
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temperatures  p.50°K  (Plate  A),  170°K  (B)  and  120°K  (C)]J  and  the  hulk 
lattice  parameter  of  4,9505  A (at  25°C),  are  4.931.  4.933  and  4.927  A, 
respectively!  also  being  in  reasonable  agreement  with  the  observed  values. 

In  contrast  with  the  above  phenomenological  explanations,  an 
opposit  effect  due  to  the  surface  atoms  may  be  predicted  on  the  lattice 
parameter.  According  to  the  calculations, the  interplanar  spac- 
ings  adjacent  to  a planar  free  surface  are  larger  than  equivalent  bulk 
spacings,  and  the  degree  of  interpleinar  expansion  diminishes  rapidly 
with  increasing  depth  below  the  surface.  For  instance,  in  the  case  of 
Pb,  the  percentage  of  increase  in  spacing  normal  to  the  surface  plane 
between  the  first  and  second  layers  for  (001) , (110)  and  (111)  surfaces 
is  calculated  to  be  5.5,  4.1  and  1.6  respectively. A LEED  ex- 

35) 

periment  for  nickel  (OOl)  surfaces  supports  such  a dilation  of  surfaces. 

For  micro  dust era,  however,  such  a surface  expansion  as  expected  in 

the  surface  of  bulk  materials  may  not  be  applied,  since  in  order  to 

stabilize  a microcluster  system  with  a large  portion  of  surface  atoms, 

the  cluster  structure  could  be  no  longer  simply  a micro-replica  of  the 

4 5 36) 

bulk  structure,  as  indeed  demonstrated  by  model  calculations,  ' 
where  an  icosahedral  structure  is  stabilized  relative  to  tlie  f.c.c. 
structure. 

Another  effect  on  the  change  in  lattice  parameters  may  arise  from 

crystal  defects.  Distortion  around  point  imperfection  in  simple  crystals 

has  been  calculated  for  a vacancy  and  for  an  interstitial  atomi  the 

37) 

shrinkage  of  the  crystal  volume  due  to  a vacancy  is  estimated,'"^  and 
for  a f.c.o.  crystal  the  first  nearest  neighbor  distances  around  a 
vacancy  decrease  but  the  second  nearest  nei^bor  distances  increase, 
and  vioe  versa  around  an  interstitial  atom.^^^ 
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Oonsldering  the  above-mentioned  structure  change  and/or  possible 
crystal-defect  effects,  we  have  to  regard  the  experimentally  determined 
lattice  parameter  for  micro clusters,  even  when  the  indexing  of 
Bragg  reflections  fits  the  case  of  the  perfect  crystal,  as  not  only 
a convoluted  q,uantity  averaged  over  a size  distribution^ but  also  an 
"effective"  quantity.  The  physical  meaning  of  the  lattice  parameter 
becomes  obscure.  The  closest  approach  distance  (r.) , calculated  from 
the  experimental  lattice  parameter  (a^)  (for  the  perfect  f.c.c.  crystal 
r^,  =ag/2^),  is  ambiguous,  and  should  no  longer  be  a simple  function 
of  Sg.  A more  proper  physical  quantity  for  r^  is  defined  by  the  first 
peak  of  the  radial  distribution  function,  although  so  far  it  has  not 
been  done  experimentally.  Model  calculations  by  Briant  and  Burton, 
using  the  radial  distribution  functions,  show  the  closest  approach 
distance  increases  with  decreasing  cluster  sizei  opposit  to  the  case 
of  the  experimentally  reduced  r^  from  the  "effective"  lattice  parameter. 

Therefore,  care  should  be  taken  for  an  interpretation  of  the 
effective  lattice  parameter.  It  is  merely  a qualitative  measure  whether 
the  microoluster  possesses  the  characteristic  of  the  bulk  crystal 
structure.  For  more  detailed  information  concerning  the  cluster  structun 
we  must  rely  on  the  analysis  of  the  Bragg-peak  intensity,  as  demonstrated 
in  the  present  study  and  further  discussed  below.  A more  elaborate 
analysis  would  be  a simulation  of  the  whole  diffraction  intensity 
contour,  using  a cluster  structure  model  based  on  molecular  dynamics 
calculations,  at  the  high  cost  of  computation. 

The  observed  periodic  deviation  of  the  Bragg-peak  intensities 
from  the  pseudo-kinematical  expression  have  been  interpreted  in  terms 
of  a partially  amorphous  structure  of  clusters.  This  may  be  justified 
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by  considering  possible  vestiges  of  liquid  cluster  structure  from  the  early 
stages  of  formation,  followed  by  a relatively  quick  adiabatic  cooling,  and  r 
large  portion  of  surface  atoms  which  may  have  many  crystal  defects 
(irregular  atomic  configurations).  Puthermore,  according  to  theoretical 
predictions,^’^^^  clusters  possess  liquid- like  surface  regions  but 
solid  central  regions.  When  the  liquid  phase  is  quenched,  a disordered 
solid  having  a structxire  similar  to  that  of  liquid  is  produced. 

However,  other  explanations  might  be  possible  for  the  diffuse  scattering. 
In  the  present  analysis,  multiple  scattering  corrections  were  made 
by  the  simple  two-beam  approximation,  and  the  thermal  vibration  effects 
were  treated  within  the  framework  of  klnematlcal  theory.  As  mentioned 
before,  for  the  higher  order  reflections  the  present  two-beam  theory 
was  certainly  insufficient.  In  multiple  scattering  situations,  the 
temperature  effects  are  more  complicated  than  the  klnematlcal  factor 
e“^i  correlation  between  displacements  on  neighboring  atoms  can  enter 
into  the  theory.  In  addition,  for  surface  layers,  enhanced  amplitudes 
of  vibrations  (of  the  order  of  50JS  greater  than  in  the  bulk)  and  their 
anisotropic  nature  have  been  predicted  In  bulk  materials.  ' These 
effects  may  also  introduce  a systematic  error  In  the  analysis.  Further 
effects  due  to  inherent  scattering  processes,  such  as  Inelastic-elastic 
scatterings  (Klkuchi  pattern),  polarization,  and  reflection  of  electron 
waves  by  the  inner  potential  of  the  crystal,  were  not  considered  in 
the  present  analysis.  On  the  experimental  side,  a size  distribution 
of  clusters  was  not  known.  The  deconvolution  of  the  observed  intensities 
by  the  size  distribution  function  folded  with  the  selective  sampling 
effect^^^  of  electron  diffraction  was  not  carried  out. 

In  spite  of  many  unsolved  problems,  the  present  interpretation 
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seems  reasonable  (Ockham's  razor),  and  lends  experimental  support  for 

the  theoretically  predicted  "liquid- like"  solid  (coexistence  of  regular 

lattices  and  random  configurations)  for  mi croc lust era?  ' If  the  present 

finding  is  inherent  to  a small  size  of  clusters,  we  could  expect  such 

a characteristic  fluctuation  of  Bragg  intensities  as  shown  in  Pig.  5 

for  other  systems.  It  is  indeed  found  in  the  case  of  Ar  clusters  with 

0 

500~'1000  atoms  (sizei  55/^42  A diameter,  similar  to  the  present 
case  (Plate  C)),  studied  by  Purges  et  al.^^  They  observed  the  anomaly 
in  the  peaks  (111),  (400)  and  (642),  suggesting  noncrystalline  effects 
of  Ar  clusters.  Using  their  observed  peak  heights  and  their  kinematical 
intensities  (given  in  their  Table  III),  the  corresponding  analysis  to 
Plate  C (Pig.  5)  has  been  made  and  is  shown  in  Pig.  6|  compare  this 
with  Pig.  5.  A striking  similarity  can  be  observed  for  both  "f.c.c." 
clusters.  By  use  of  eq.  (14)  with  <r  * 3.4  A,  t|^  = 0.45,  and  2M  = 0.021s^, 
the  characteristic  fluctuations  observed  in  the  Ar  clusters  can  be  well 
reproduced,  as  shown  in  the  upper  curve  of  Pig.  6. 
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Table  1.  Experimental  Conditions  for  Pb  Cluster  Beams. 


Plate 

Tn(°C)®^ 

Ta(°C)'*^ 

Pa(torr)°^ 

Exp.(sec)®^ 

A 

993 

79 

0.83 

3.6 

30 

B 

1040 

150 

0.67 

2.8 

20 

C 

932 

103 

0.55 

2.8 

20 

a)  The  tmiperature  of  the  tungsten  filament,  measured  by  an  optical 
pyrometer!  uncertainties  of  c.a.  ± 100°C.  The  electric  power 
required  for  these  twnperatures  are  12^14  Amp  and  -^1  V, 

b)  The  temperature  of  the  mixing  region,  measured  by  a Chromel-Alumel 
thermocouple!  see  Pig.  1. 

c)  Argon  carrier  gas  pressure  in  the  mixing  region. 

d)  The  distance  between  the  heater  (H)  and  the  first  nozzle  aperture 
(N^) I see  Pig.  1. 

s)  Photographic  exposure  time. 


i 


I 


20 


P 

V 

>, 

P 

i 


» 

I 

i 
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Table  2.  The  observed  lattice  parameters,  Bragg-peak  intensities, 
average  cluster  sizes  and  cluster  temperatures  for  Plates 
A,  B and  C. 


a)  Observed  peak  intensities,  normalized  to  100  at  the  (ill)  peek. 

b)  Average  of  all  observed  lattice  parameters!  uncertainties  are  2.3 
times  their  standard  deviations. 

c)  Average  cluster  size  in  diameter,  estimated  from  the  line-vidth 
analysis!  see  text. 

d)  Cluster  tmperatures,  estimated  from  the  Debye-Waller- factor  analysis! 

see  text. 
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Figure  1 
Figure  2 


Figure  5 


Figure  4 


Figure  5 


FIGURE  CAPTIONS 

Schematic  diagram  (not  to  scale)  of  the  nozzle  for  metal 
cluster  beams  I see  text. 

Photographic  density  (O.D.)  of  the  diffraction  patterns  of 
Pb  clusters  for  Plates  A,  B and  C,  plotted  as  a function  of 
the  scattering  variable  si  Miller  indices  are  given  in  A. 
Logarithmic  plots  of  the  ratio  of  the  observed  Bragg-peak 
intensity  1^^^  (Plate  A)  to  the  calculated  intensity  I^aic* 
plotted  as  a function  of  s^  (Debye-Waller-f actor  analysis). 
Icaic  ^dyn  text)  was  calculated  using  partial 

wave  (PW)  and  Born  atomic  scattering  factors  with  various 
size  parameter  Dt  (1)  D ^ 90 A (PW) » (2)  D — 80A  (PW)  i (3) 

D = 60A  (Born)  I (4)  D ~ 60A  (PW)  i and  (5)  D « 0 (PW) , equiv- 
alent to  Igalc  ~ ^kin  ®^*  ^®**  (kinematical  analy- 

sis) • 

Plots  similar  to  Figure  3 (Debye-Waller- factor  analysis) 
for  Plates  A,  B and  C.  Only  the  best  straight  line  for 
each  Plate  is  showni  for  A the  same  as  (2)  in  Figure  3,  for 
B D * 60A  (PW) , and  for  CD*  50A  (PW) , respectively. 

The  fluctuations  R(s)  of  Bragg-peak  intensities  (Pb  clusters  1 
Plate  C) . The  lower  curve  is  the  observed  R(s),  defined  by 
eq.  (8),  and  the  upper  curve  is  the  calculated  R(s)  by  eq.  (14), 
based  on  the  "liquid  model"  described  in  text.  The  calculat- 
ed points  are  normalized  to  1.28  (the  observed  value)  at  (111). 
The  dotted  lines  are  merely  drawn  to  accentuate  the  charac- 
teristics of  fluctuations. 
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Figure  6.  The  fluctuations  R(s)  of  Bragg-peak  intensities  for  Ar 

clusters  I plots  similar  to  Figure  6.  The  observed  points 
are  taken  from  the  Table  III  of  Ref*  li  see  text.  The  calcu- 
lated points  are  normarited  to  1.5  (the  observed  value) 
at  (111). 


